Dynamical scaling of the critical velocity for the onset of 
turbulence in oscillatory superflows 



R. Hanninen W. Schoepe 



March 25, 2011 

Abstract The critical velocity iv for the onset of turbulence in oscillatory flows of super- 
fluid helium is known to depend on the oscillation frequency CO, namely vv ~ \/kco where 
K is the circulation quantum. Only the numerical prefactor may have some geometry de- 
pendence. This universal behaviour was described earlier qualitatively either by employing 
the superfluid Reynolds number or by extending known dc vortex dynamics to ac flow. In 
our present work we emphasize that vv(a)) °<: ^/a can also be derived rigorously by means 
of dynamical scaling of equations of vortex dynamics as pointed out by Kotsubo and Swift 
already two decades ago. 
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1 Introduction 

Dynamical scaling of an equation of motion is a simple but powerful tool for deriving rela- 
tions between various parameters of that system. For example, in a classical viscous liquid 
the Navier-Stokes equation 

d,\+{\-V)\ = -^Vp + vV\ (1) 

is invariant when lengths are scaled by / = A /*, velocities by v = v*/A, and times by ? = 
[1]. This scaling also leaves the Reynolds number unchanged 

Iv l*v* 
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which is the famous similarity principle of fluid dynamics. 

In the following we will summarize what is known about dynamical scaling of the motion 
of vortex lines in superfluid helium. 



2 Scaling of vortex dynamics 

In his seminal paper of 1988 Schwarz [2] has scaled his equation of motion of a vortex fila- 
ment s = s(<§,f) 

s = jSs' X s" + V.5 + as' X (v,„ - jSs' x s") 

-a's'x[s'x(v,„-jSs'xs")], (3) 

where the primes denote derivatives with respect to the arc length <§, a and a' are the 
parameters of mutual friction, and jS = {K/An)\n{cR/aQ), with c ~ 1, is an effective 
radius of curvature of the filament, and uq the core radius. The counterflow velocity, v,„ = 
— Vj, is simply the difference between normal and superfluid velocities. After introducing 
reduced time and velocity scales to = fit and vo = I'/jS the resulting equation was shown to 
be invariant under the same scaling as discussed above. 

Even if the above scaling was derived by using a so called local induction approximation 
(LIA) the same applies also when using a full Biot-Savart equation (BSE). The inclusion of 
the non-local term, which itself satisfies the velocity scaling, only slightly modifies the local 
term. Since, as will be noted below, the corrections due to the logarithmic term are small and 
very difficult to observe in experiments, we ignore here the small changes in the logarithmic 
term that are due to differences between LIA and BSE. 

It should be mentioned that from Eq.(3) Schwarz calculated a rate equation of the change 
of the vortex density dL/dto due to a growth term L^^^ and a loss term Lr that is scale 
invariant and that corresponds to the well known Vinen equation [3], see below. 



3 Scaling of oscillating superflows 

In ac flow experiments the oscillation frequency o is a scaling variable. This has been rec- 
ognized some time ago by Kotsubo and Swift when analyzing their experiments on vortex 
turbulence generated by high-amplitude second sound in helium-4 [4,5]. In the following, 
we apply their scaling results to our experiments on ac flows generated by an oscillating 
sphere. From the work of Schwarz [2] they find that the time scale 1 /co transforms as 

fflo = A--(Bo*, (4) 

where coq = co/P and cOg = ro*/jS* with jS* = (k /47i) In {cR* /ao). From v^o = ^^^v*o they 
obtain by eliminating A as a final result 

vv/v: = (/3//3*)i/2(a)/co*)'/2. (5) 

This is how the square-root dependence of Vc{(x)) is derived. There is, however, the logarith- 
mic con'ection by (jS /li*y^^ due to the scaling R = ?lR*. We can estimate the logarithmic 
corrections from the definition of /5 and jS* and using Eq.(4) by solving for jS/jS* [5] 

= (6) 

In (cR'/ao) 
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If we assume that R and R* are approximately given by the intervortex spacing as suggested 
by Schwarz [6] we estimate as an order of magnitude In (cR* /ao) w 10. Inserting this we find 
(j6/j6*)'/^ = 1.04 for our complete frequency range a* /(O w 7. This is a small correction 
that could not be resolved in our experiment, see Fig. 1 [7]. 
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Fig. 1 (From Ref.[7]) Critical velocity for the onset of turbulence as a function of the frequency of 2 oscillat- 
ing spheres. Blue squares: radius 124/im; red dots: radius 100/Jm; green soUd Une: iv = 2.85 \J Klnf (m/s); 
temperature is below 1 K. 



It is true that in our experiments the sphere radius is kept fixed and therefore the square- 
root dependence of is not an automatic result from scaling. But as noted by Kotsubo and 
Swift, the scaling suggests three possible scenarios [5]. One is that Vc depends only on the 
frequence via Eq.(5), as appears in our experiments. Second option is that it depends only 
on the geometrical dimension D. This seems to be the case for dc flows as summarized, 
e.g., in Refs.[8,9]. The last option is that could depend on both O) and D. In this case the 
square-root dependence would show up only in a particular experiment where the frequency 
and dimensions are both scaled properly. 

While the frequency dependence of Vc can be found from scaling no information on 
the absolute values or the temperature dependence can be inferred. As pointed out earlier 
[10] the prefactor 2.85 of our experimental result on oscillating spheres Vc = 2.85 \/ JCO) (see 
Fig. 1) may have some geometry dependence because the superfluid velocity field in general 
varies over the surface of the oscillating body. Hence, this prefactor can only be determined 
from experiment. In fact, the vibrating wire data of the Osaka group give a slightly smaller 
value [11]. On the other hand, iv from the second-sound data of Kotsubo and Swift is about 
twice of what we find from our spheres at the same frequency and temperature. 

Another difference concerns the temperature dependence of vv. Kotsubo and Swift find 
that their Vc decreases with rising temperature from 1 .2 K to 2.0 K, while in our case we 
observe a 10% increase from 1.0 K to 1.9 K that follows from the Vinen equation, see 
below. Possibly, these differences may be due to the very different mechanisms by which 
the oscillating flows are excited. 
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4 Conclusions 

In our earlier work [7, 10, 12] we have presented two different ways of deriving the fre- 
quency dependence of A qualitative but very general argument was based on the super- 
fluid Reynolds number Re^ = Iv/k ^ 1 by identifying the oscillation amplitude v/co as the 
characteristic length scale /, which immediately implies the square-root behaviour of Vc(co). 

Secondly, and in more detail, we started from the Vinen equation [3], which was estab- 
lished for dc flow and recently modified by Kopnin [13] 

dL/dt = b{vL^I^ -kL-), (7) 

where b=A{\ — a') —Ba and A, 5 ~ 1. (Equation (7) does not contain the logarithmic term 
jS and is invariant under the same scaling as given in the Introduction.) Applying the Vinen 
equation to ac flow we have postulated that the vortex tangle, in order to follow the amplitude 
V of the ac flow field, must have a relaxation time T = iK/bv [13] that is shorter than the 
period of the velocity field, i.e., a)T < 1. This gives Vc(a)) ^ yJlKCoJb. Below 1 K mutual 
friction can be neglected, hence b ~ \. Towards higher temperatures b gradually decreases 
and hence increases. This is in agreement with our data [12]. Clearly, both methods have 
a qualitative character. But to our knowledge, a theory of oscillatory superflow, in particular 
an extension of the Vinen equation to ac flow is not yet available. 

In contrast, our present approach is based on dynamical scaling that follows from sym- 
metry properties of the equations of vortex dynamics. Equation (5) is a rigorous result. The 
logarithmic terms contain the scaling of the length which in our experiments, however, can 
be neglected. Therefore, the square-root dependence of Ve(a)) is well established. 
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